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I. INTRODUCTION 

The phenomenological appeal of planar QCD is well known [Q . It is believed that a limit in which the number of 
colors iV c goes to infinity while the coupling goes to zero with fixed g 2 N c produces the simpler planar theory while 
still preserving essential perturbative and non-perturbative features of the physical N c = 3 case |^|, However, an 
analytical description of the planar limit has remained elusive. This has left lattice gauge theory as the prevailing 
approach for the conceptual and quantitative analysis of strong coupling in QCD. 

Recent numerical investigations with increasing N c |^] indicate that some N c — 3 physical quantities arc not 
far from their N c = 00 limits. This suggests that numerical results from planar QCD might have phenomenological 
value. However, in numerical lattice calculations, it is easier to deal with the N c — 3 case than to approach the planar 
limit by simply increasing N c . An approach that arrives at the simpler planar theory more directly is indicated. The 
Eguchi-Kawai |(| reduction offers such a shortcut. It has the potential to require much less computation because 
the planar limit of some observables is obtained from a matrix model with only the d x i\T? degrees of freedom of d 
unitary matrices. Since the dimension d enters linearly rather than exponentially, it seems that four dimensions has 
only twice the computational cost of two dimensions. 

Using the Eguchi-Kawai reduction to approach the planar limit of QCD by first going to infinite N c and then to 
zero lattice spacing might be beneficial because the two limits may not commute [0, |g] . Also the commutativity of the 
infinite volume and infinite N c limits might be questioned because gauge theories are massless, and the construction 
of H docs not apply. We begin with the assumption that meson-meson correlation functions at non-zero momenta 
are "infrared safe" . If our calculations in two dimensions are finite and agree with analytical results, then we can be 
confident about calculations in four dimensions, where infrared problems are less serious. 

In the Eguchi-Kawai model the entire d-dimensional lattice is reduced to a point. Nevertheless the matrix model 
incorporates the lattice structure. This is best understood in momentum space ]jfj| , where the d sets of eigenvalues 
of the link matrices populate the momentum d-torus. To ensure that the torus is indeed covered, one needs to alter 
the Eguchi-Kawai proposal. The first suggestion for how to achieve this was to freeze the eigenvalue variables and 
remove them from the integration over gauge fields. The fluctuations of the remaining gauge degrees of freedom are 
not allowed to back-react on the eigenvalue distribution. The expectation value of gauge invariant observables is 
subsequently averaged over different sets of angles. This produced the quenched Eguchi-Kawai model jl(|. Related 
work on quenched reduced models appeared in 111], |12l O]. A more elegant method, which involves a change in 
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the action by a "twist" [|14|, was subsequently discovered. This clever approach has also led to the discovery of 
non-commutative gauge field theories [pL 5| starting from the lattice. 

Nevertheless our proposal is in the quenched version jlO| because the introduction of fermions in the twisted case is 
somewhat restricted. The addition of fermions to reduced models was first carried out in a continuum version in [jn]. 
A somewhat similar procedure on the lattice was introduced in In the latter version, the number of flavors 

Nf was a free parameter and could be taken to infinity. This produces a reduced model leading to the Veneziano limit 
p8[ rather than to the 't Hooft limit jlj . Following up on the idea of [[l6] , a twisted version of the Nf = N c case was 
presented in [ist| . Since we want to have a variable number of flavors, we opt for the quenched version. In addition, 
the quenched version is computationally simpler than the twisted one even for Nf = N c because the Dirac operator 
is block-diagonal with Nf blocks each of size N c x N c . In the twisted case, one has a full N c Nf x N c Nf matrix. 

Our proposal specifically targets mesonic observables. To get glueball masses from a reduced model, we would 
need to keep one direction unreduced as described in p0| . To investigate finite temperature, we would again need a 
Hamiltonian-like formulation pl| . In addition to the calculation of meson masses, our approach allows for the study 
of fermion quenching effects (valence approximation effects) and finite fermion number density at zero temperature. 
Also the pathologies of the valence approximation should disappear in the limit of large N c with fixed Nf. The 
fermionic bilinears which define the meson fields are the building blocks of the four-fermi operators crucial to the 
study of weak decays. We envisage taking the ^- approach to non-leptonic weak decays j2^] to the lattice where it 
could be compared to more direct approaches and where some of its phcnomenological input could be replaced by 
numerical data obtained from "first principles" . 

There have been both practical and theoretical developments that motivate a return to numerical work on planar 
QCD. Theoretical progress has led to the new lattice fermions p3|, which might be able to solve some problems 
encountered about twenty years ago |TT) . The problems had to do with the meaning of topological charge density 
in the reduced model, chiral anomalies, and the ^-parameter dependence of the vacuum action density in the planar 
limit p4fl . Equally important are developments in computer hardware. Since the initial numerical experiments on 
reduced models for planar QCD, computers have become 10 4 times more powerful for fixed power consumption. We 
suggest that reduced model simulations are a worthy job for weakly (and cheaply) linked PC clusters. Moreover at 
infinite N c , the Wilson Dirac operator in the argument of the sign function has no small eigenvalues when the pure 
gauge action of the reduced model is of the usual single plaquette type, and the rescaled lattice gauge coupling is large 
enough (but not infinite). This is unlike the situation in ordinary QCD p5[ . Thus the numerical implementation of 
the overlap Dirac operator will be simpler than in ordinary QCD. 

In the next section, we define reduced models starting from ordinary formulations of lattice gauge theories. We 
show that the topological charge remains a nontrivial object even after reduction and that reduced models have, in 
some sense, more topological numbers than their unreduced counterparts. We briefly discuss chiral gauge theories. 
The meson-meson correlators, which are the most practical aspect of our proposal, are also introduced. Section three 
is focused on the special case of two dimensions, in which we test the numerical methods. Although planar QCD2 
has been solved exactly in the weak coupling phase, additional analytical and numerical work is needed to extract 
the results that can be compared to Monte Carlo data. Although there is a vast literature on QCD2, in the interest 
of clarity, we have included a rather detailed account in the Appendix. In the summary section, we describe some of 
our ideas for future developments. We believe that there are many options for future research and that the potential 
for truly new and interesting results is quite high. 

II. GENERAL FRAMEWORK 
A. Original lattice model and its topological charge 



Before reduction, the lattice systems of interest are defined on a d-dimensional hypercubic lattice with a total V 
sites on a d-torus of equal sides. Sites are labeled by x and positive directions by fi = 1,2, ... ,d. There are SU(N C ) 
gauge fields U fl (x) on the links, and fermions interact with this gauge background. 

Since we are going to deal later with fermions in the fundamental representation, it is best to introduce the fields 
now, even though we are first focusing on the pure gauge sector. The fermion field ijj is a vector in V s ® V c <8> Vi , the 
tensor product of the spin, color, and lattice vector spaces, respectively. Its components are ip l a {x) with spin index 
a, color index i and lattice position x. There are several operators acting on fermions: The Euclidean Dirac 7 M 's act 
only on spinorial indices, 7^ : V s — > V s . The directional parallel transporters T M act on the site index and the group 
index (T M : V c <£> VJ — > V c <g> Vi ) as follows: 

T^)(x) = U^x)^(x + fl). (1) 
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For Ufj,(x) = 1, the T M become commuting shift operators T°. For V = L d each T° is unitary and has eigenvalues 
e 2 f ife , k — 0, 1, 2, . . . , L — 1. Hence for even d, T M is not only unitary but also has unit determinant (detT° — 

Gauge transformations are characterized by a collection of g(x) € SU(N C ) acting on ip pointwise, and only on the 
group indices. The action is represented by a unitary operator G(g) : V c ® Vi — > 14 ® Vi with (G(g)^)(:c) = g(x)ip{x). 
The T M operators are gauge covariant 

G(g)T ll (U)Gt( g )=T ll (U*) (2) 

with 

C^( a! )= ff (a:)t^(x) ff t( a; + /i). (3) 

The variables U^x) are distributed according to a probability density that is invariant under U — > U 9 for any g. 
For the commutators [T^, T u ], we have 

[T^T^T^] = (1 - P^)t(l - P pi/ ) = 2 - iV - Pt„ (4) 

with the unitary given by 

Pfiv = TjjT^TJT^. (5) 

The operators P^„ are site diagonal with entries that are the parallel transporters around plaquettes: 

(P^)(x) = U^ v {x)^{x) (6) 
U MV (x) = Ul{x - 0)U^(x — v- jl)U v {x -i>- p,)Up(x - jl). (7) 

Ufj, v {x) is associated with the elementary loop starting at site x, going first in the negative v direction, then in the 
negative /x direction, and coming back round the plaquette. 
The Wilson pure gauge action is given by 

S g = [3 Tr(P^ + Pt „) = Const - £ Tr[T M , T^[T^ T v ] . (8) 

In the continuum limit, configurations with all the U^ v (x) close to the unit matrix are preferred. 

All our operators are finite dimensional matrices. The norm of A is \\A\\ with ||A|| 2 the largest eigenvalue of A* A. 
We have 

\\[T il ,T v ]\\ = \\i-p liV \\. (9) 

The gauge invariant constraint 

\\[TM\\<e^ (10) 

for all /i > v does not change the continuum limit. It is equivalent to 

\\l-U^{x)\\ <e M „ (11) 

for every site x. The are small fixed nonzero numbers, independent of (3 and of N c . The probability distribution 
of the link variables is e s with integration ranges restricted by ([l0|) . 

The lattice version of the massive continuum Dirac operator D(m): V s ® V c <8> Vi — ► V s ® V c ® Vi is an element in 
the algebra generated by T M , T^, 7^ and thus is gauge covariant. The Wilson Dirac operator Dw(m) is the sparsest 
possible analogue of the continuum massive Dirac operator which obeys hypercubic symmetry. Fixing the r-parameter 
to its preferred value (r — 1), 

£V = m + £(l-X M ) XlX, = l X ,= l —^T,+ l -±^Tl. (12) 
One can then prove a general bound |2q| : 



m in 



> 



l-(2 + V2)^ £ 



-|l + m|. (13) 
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Here A m i„(M) is the smallest eigenvalue of the matrix M. This bound has content only for m close to —1: 

i 

1 2 

I in •_ l-(2 + V2)J2 € ^ ■ 



(14) 



This range is contained in the open interval —2 < m < 0. 
For even d, we define the hermitian operator 

H w (m) = j d+1 D w (m) for which H^(m) — D\ v (m)D w {m). (15) 

It follows from © and (H) that H w has no zero eigenvalues so that its sign function e = sign(_ffw) is well-defined [fj7| . 
The topological charge of a gauge field configuration is given by 

Q = \Tre (16) 

The value of Q cannot be changed by a smooth deformation that respects the bounds ([!(]) . Close to continuum, Q 
equals the continuum topological charge j2^]. The trace in (|l6|) includes a sum over lattice sites. Leaving this sum 
out produces a local functional of the lattice gauge background which is a lattice version of the topological charge 
density trFF(x). For a review, seeJ29j. 

It is important that the proof [^6| of the bound makes no use of the special structure of the matrices 7), in (|l|) 
except their unitarity. Thus any unitary 2), with commutators close to zero ( jioj ) can be used in (|l2] ) and (|l^) to 
define homotopy classes for the gauge backgrounds. Of course, the "local" topological density has no meaning if the 
T^'s are that general. 

Although we spoke about fermions and used gamma matrices, they were auxiliary concepts. We only dealt with 
the pure gauge sector. The gamma matrices played a role in defining the topological number, but the latter is just a 
function of the gauge field background. It will be important in the following discussion that this definition of Q can 
be non-zero on a finite lattice. One need not to go to an infinite lattice to see non-zero topology. 



B. Reduced model and reduced topological charge 

The reduced model for lattice gauge fields || is defined on a toroidal lattice with a single site. The matrices T M 
become unitary matrices with no restrictions on their structure since they are equal to (cf. (0)) on a single site 
lattice. A connection between models differing only in the number of sites V can be established by taking N c to infinity 
with (3/N c fixed. The leading behavior of traces of Wilson loops written in terms of the T M is then independent of V. 
Actually this is only true for j3jN c small enough. It fails for large j3/N c and thus also for the continuum limit |L0| . 
The key to fixing the model is to understand that the V independence is caused by momentum space getting absorbed 
into the space of eigenvalues of the operators T M . Consequently, the quenched reduced lattice model is defined by 
freezing the eigenvalues of T M to a set uniformly distributed over the unit circle for each /i ]ic| . 

T M = VpDjQ, = diag(e< , e< , • • • , e<° ) (17) 

The angles O 1 ^ are randomly distributed in the interval [— 7r, 7t]. The ordering of the angles for each direction \x does 
not matter since one can go from one ordering to another by appropriately changing the integration variables . In 
other words, only the sets of eigenvalues are frozen, not the ordering within each set. The matrices are unitary 
and otherwise restricted only by ([To|). 

The pure gauge action is invariant under — > V^D'^ where the D' are arbitrary diagonal unitary matrices. Thus 
the Vfj, should be thought of as taking values in the coset space U (N c ) / (U (1)) N " . The most important symmetry of 
the pure gauge action is gauge invariance which now acts by conjugation T M — > flT^il^ . One can think of f2 as the 
replacement for G(g) in (||). The gauge group is SU(N C )/Z(N C ). 

The absolute minima of the action are given by T M 's that are simultaneously diagonal. These minima are changed 
one into another by certain sets of matrices. In the double line notation [[[o| [ll[ [l2[ |l3) of planar graphs, the 
crystal d-momentum going through a line is given, component by component, by the difference between the angular 
phases of the two eigenvalues of T M associated with each of the color indices attached to each component of the double 
line. Different minima give different momenta, and the sum over minima together with the sum over color indices 
reproduces the ordinary planar lattice Feynman diagrams. As the typical spacing between a O 1 ^ and the closest 6^ is 
of order l/N Cl it follows that the momenta values are in effective correspondence with those in a finite volume of size 
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V = N c . On the other hand, focusing on a single minimum might suggest V — N c . In the twisted case, one typically 

d 

gets V — N c 2 . In numerical simulations, it is important to ensure that all minima are properly sampled. 

In the continuum version of the reduced model pj] j, the basic variables are traceless hermitian matrices with 
restrictions on their eigenvalues. One could try a "geometric" definition of topological charge, but it does not work 
well O] . The natural expression corresponding to the non-abelian field strength is F M „ = [A^ , A v \ . So long as the 
A„ are finite matrices, 



TrFF ex euv P <rTrAaA v A p A a = 



(18) 



because of the cyclicity of the trace and the antisymmetry of the epsilon symbol. To get non-zero topological charge, 
one needs to work with infinite A^ matrices, and a numerical approach seems impossible. The topological charge of 
the lattice reduced model is defined by ( |l6|) using unitary X^'s that are restricted only by ([n]). As already mentioned, 
the definition of topological charge and the restriction on the commutators (10) work exactly the same way as before. 

This extra constraint (10) on the commutators of T M is not necessary at infinite N c . For d = 2, when we get close 
enough to the continuum limit and the lattice gauge coupling b 2 — j3/N c is large enough, we are far below the Gross 
Witten phase transition |30| , and the plaquette eigenvalue distribution is contained in a small neighborhood of unity. 
This is evident in Fig. I] where we have plotted the expectation value of the maximum angle of the untraced plaquette 
in units of 7r for two different values of iV c in the quenched model. This is also true in four dimensions where again 
there is a phase transition as seen in the Fig. [l]. A similar result was also found in the four dimensional twisted 
Eguchi-Kawai model ]3l| . The numerical evidence strongly indicates that the extra condition restricting the matrices 
is satisfied as long as we are close enough to continuum even if we use only the simplest single plaquette Wilson action. 

A complete study of topological charge is beyond the scope of this paper and our discussion will be brief. Consider 
a two dimensional gauge field background in which the T M 's represent a U(l) background that is a U(l) instanton on 



a two dimensional torus. 



The T M have the structure of a 



'shift" and a "clock" matrix l3z 



Ti = X 









V 1 





1 

0/ 



(19) 
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Here ( = e~ , and we take N = N c . This abelian background on a lattice of area N c can be viewed as a non-abelian 
SU(N C ) background in the reduced model. As far as the spectra of the various lattice Dirac operators go, how we 
look at the background is irrelevant. Hence, the reduced model has Q — 1 for this configuration. Moreover calculating 
the gauge field action, we easily can check that it stays finite in the 't Hooft large N c limit. Similarly, one can easily 
check that ([!(]) is satisfied for N c beyond some finite value. 

One expects no topological effects in two dimensional planar QCD. This is not in contradiction with our observation 
because the continuum limit is obtained when the fixed parameter of the 't Hooft limit f3/N c is taken to infinity, and 
these backgrounds do not contribute. We learn that lattice reduction may perform a "miracle" in the sense that one 
may find topological effects that were not present in the original continuum model. Whether or not these effects are 
relevant in the continuum limit is a dynamical question. 

What happens in four dimensions? We have seen that what looks like an abelian background before reduction can 
appear as a nonabelian background in the reduced model. (Actually, the results of |3^] say that for any finite N c > 3 
on a continuous torus, any allowed value of Q can be attained by an abelian background. We are not making use of 
that result here.) In four dimensions, one can obtain a configuration that is nontrivial topologically, for example, by 
choosing an abelian instanton in the 1 — 2 plane and another abelian instanton in the 3 — 4 plane. These instantons 
would lead to shift-clock pairs as above. We could take N c — N 2 and write the T^'s as direct products of N x N 
matrices: Ti = X ® 1, T 2 = Y 1, T 3 = 1 ® X and T 4 = 1<E>Y. All these matrices are in SU(N 2 ). The topological 
charge of this background is Q = 1. The spectra of the T^'s are highly degenerate in this example. To produce truly 
acceptable configurations in a quenched reduced model, we would need to find deformations of these matrices that 
lift the degeneracy while preserving Q = 1 . 

The action of this configuration diverges linearly with N c , just like the action of a continuum instanton. But we see 
that there are configurations with non-zero Q, and the question of their survival at infinite iV c is similar to the familiar 
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FIG. 1: Plot of the maximum plaquette angle in units of it as a function oib 2 = The data is at N c = 31 and iV c = 37 for 
the two dimensional model and is at N c — 17 and N c = 19 for the four dimensional model. 
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question in the continuum |34[]. This indicates that in four dimensions, the reduced lattice model could provide a 
framework in which one can at least pose questions about topology at infinite N c . This is interesting because the 
results of Teper 0] support a finite nonvanishing limit to the topological susceptibility in the lattice planar limit. We 
do not know if there are configurations that produce Q — 1 and have an action that stays finite as N c goes to infinity. 

We now make an observation that might be of relevance to the issue of ^-dependence in the planar limit. The set of 
backgrounds satisfying (|l^) falls apart into non-empty disconnected components beyond what is needed for distinct 
topology. Even at fixed Q there are several components. These different components make separate contributions 
to the partition function, and these contributions are not obviously suppressed at infinite N c . It could be that the 
particular connected component that dominates at infinite N c depends discontinuously on parameters, and therefore 
the action density need no longer be differcntiable with respect to its parameters. In particular, this could induce a 
nonanalyticity in 9/N c , and this might be a way of reconciling the naive scaling for the action density (~ N 2 f(^-)) 
with 2tt periodicity in 9 [f24| . 

To see that the set of gauge fields one integrates over splits into disconnected components, consider again the 
above two dimensional example. Complete the pair T\,2 by a T3 and a T4 of a perturbative nature so that Q = 0. 
The shift-clock pair cannot be deformed smoothly to unity while obeying the constraint on the commutator because 
the shift-clock pair could be used to construct a two dimensional lattice Dirac operator which would have a two 
dimensional topological charge of unity. Thus, this set of four T^'s has Q = but is homotopically inequivalent to a 
set where all four Tj/s are perturbative and consequently also has Q = 0. In spite of the fact that we are considering 
only SU(N C ) and not U(N C ), the two dimensional Chern numbers associated with the different two dimensional tori 
included in the four dimensional torus manage to survive in the reduced model. 

Let us briefly mention anomalies associated with chiral fcrmions. The free energy of the reduced model with fermions 
included is supposed to reproduce the free energy of the original model at leading order N% and subleading order N c . 
Hence, anomalies should spoil the gauge invariance of the reduced model too. In the overlap formalism, there are 
simple lattice abelian gauge backgrounds fl35| in which one can prove that gauge invariance cannot be preserved if 
there are anomalies in the continuum. Moreover, this happens already on finite lattices. These backgrounds can be 
"borrowed" to produce nonabelian backgrounds in the reduced model with a similar effect. The overlap line bundle 
is twisted over the torus defined by variables |pi,2| < tt that enter through the family of gauge backgrounds described 
by Ti = e lpi ®Y,T 2 = e lp2 <g) Y, T 3 = X <g> 1, T 4 ' = Y ® 1 ancLas in the case of the "instanton" , N c — N 2 . Again, we 
see the potential to move beyond the obstacles identified in [jllj , but there is a degeneracy of eigenvalues we need to 
eliminate. 



C. Fermions 



In the original unreduced lattice model, massless fermions are incorporated through the massless overlap Dirac 
operator |^7j 

Do = i(l+ 7 d+ic) (21) 

with e = svgn(H\y)- The operator D Q has an implicit dependence on the fermion parallel transporter. 

The reduced version of D is obtained by replacing the original transporter with the reduced fermion parallel 
transporter. For fermion momenta — 7r < p^ < n and given in Jl^), the reduced transporter is e lp ^T^. With that 
replacement, D a becomes the reduced massless overlap Dirac operator. It depends on the fermion momentum. The 
external line propagator f36[l is 

9(p) = ^-[D- 1 -l] (22) 
Zm 

Consider the connected correlation function of the fermion bilinear -j={p(x)ip(x). In momentum space it is given 

by 

G(P) = - — {[(Trg(p)g(q)) 

gauge average J /angle average (23) 

with P = p — q. Similar expressions can be derived for the correlation functions of other bilinears. 

The way reduction works is best understood by looking at an example. Figure ^| is a diagram showing the various 
momenta in the one gluon exchange contribution to a scalar-scalar correlator. All crystal momenta come from the 
diagonal entries of the fermion transporter and flow in the directions indicated by the arrows. The 6 momenta are 
associated with color and the p, q momenta with flavor. Single lines are fermions, and the double line is the exchanged 
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FIG. 2: A planar diagram representing the one gluon exchange contribution to the scalar - scalar correlation at momentum P. 

gluon. The gluon carries momentum 0j — Oj. With the quenched 0j randomly distributed on the momenta tori, the 
summation over the color indices i and j induces the needed integrations over momenta. 

Note that we compute the meson correlators directly in momentum space. This is an important new feature with 
several potential applications. 



When the quarks are massive, the propagator g(p) is 37 



2m (1 + fi) + (1 - Mi7d+ie 

The parameter /i fixes the bare lattice quark mass via m q = 2m/j. 

The transporter e ip "T fl introduces a phase into the reduced fermion operator. This suggests that there are several 
options for defining Q, something we glossed over in the previous subsection. It is a bit strange that one seems 
to be able to define a topological charge that depends on the momentum carried by the fcrmions. Probably, for 
configurations typical to the continuum limit, Q ends up not depending on p. This is indeed the case for examples in 
section IIB as long as N is large enough and the momenta are not too large. 

It is quite easy to add Nf flavors of fermions and then go to the topological limit Nf, N c — + oo with Nf/N c kept 
fixed. In this limit, the fermionic determinant needs to be included in the gauge field integration. Similarly, one can 
also introduce fermions in representations other the fundamentals. Chiral gauge theories will have anomalies because, 
as mentioned before, the backgrounds that produce them have a reduced model representative pq|. 



III. TWO DIMENSIONAL PLANAR QCD 



As an exercise, we will calculate the scalar-scalar 

(4>(x)x(x)x(0M0)) (25) 

and the pseudoscalar-pseudoscalar 

$(i)73X(a0x(a07J»lKO)) (26) 

correlation functions in momentum space in the planar limit of two-dimensional QCD. (We use two different fields, 
tp and x to eliminate the disconnected piece.) These correlation functions are known exactly and are completely 
determined by their Feynman diagrams. Nevertheless, they are quite nontrivial. While they can be expressed as an 
infinite sum over stable mesons with mass squares that are evenly spaced asymptotically, they also have free field 
behavior at large momenta. Thus, we have the cohabitation of asymptotic freedom and total confinement that is also 
present in four-dimensional planar QCD. Moreover, with an appropriate choice of limits, chiral symmetry is broken, 
and there a Goldstone boson. If the potential infrared problems turn out to be tractable in two dimensions, we can 
be confident that there will not be any difficulties in four dimensions where infrared issues are less serious. 

In two-dimensional QCD, spontaneous chiral symmetry breaking only occurs if one takes the limit m q — > after 
the planar limit N c — > oo at fixed j3/N c . In the reduced model for any finite N c , we have exact chiral symmetry 
at m q — 0, and we can investigate the planar limit of the massless reduced model. The continuum limit is taken 
last. With m q set to zero before the large N c limit, we will not see direct signals of spontaneous symmetry breaking. 
The simplest assumption is that chirally invariant quantities will be well-described by 't Hooft's solution, but chirally 
noninvariant observables will not be. Thus the invariant combination to compare to the 't Hooft solution is 



$(aOx(aOx(OW(0)> - $(z)73X(aOx(0)73^(0)) . 



(27) 
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For massive quarks, the situation is less subtle, and one can consider the correlators separately. The planar limit of 
the reduced model is supposed to produce a UV regulated approximation to the 't Hooft solution. It is interesting to 
see in detail how this works. On the one hand, the very high mass mesons are meaningless entities at finite ultraviolet 
cutoff, while on the other hand, they are responsible for the free field short distance behavior. 

It is believed that ordinary, massive two-dimensional QCD has two phases [ p8| : a weak coupling one in which the 
gauge couping divided by the quark mass is small and a strong coupling one in which this ratio is large. In the strong 
coupling phase, there are light baryonic states, but there is no pion whose mass vanishes as y/m q in the chiral limit, 
and it is not clear how the large N c limit should be defined, 't Hooft 's planar solution is in the weak coupling phase 
and does have a pion mass that vanishes as ^Jm q . It is interesting that the numerical results obtained in j^] seem 
to correspond to the strong coupling phase and, nevertheless, also seem to admit a large N c limit with the 't Hooft 
scaling of the gauge coupling. Thus it could be that in two-dimensional QCD, there exists a strong coupling planar 
phase that is defined through the lattice and which differs from 't Hooft's solution. 

We will work in the weak coupling phase where 't Hooft's solution should be valid. If there are indeed two planar 
phases, convergence to the planar continuum limit of the weak coupling phase might be slower than anticipated. 



A. Pure gauge 



We can go to a "unitary gauge" in which T2 is fixed to the clock matrix Y, with N = N c . Denoting T\ by U the 
partition function is 

Z=f dWe s S = 4/?V|iy 2 sin 2 ^— A. ( 28 ) 

JU=WYW, W€U(N) y N c 

Terms with i = j do not contribute to S. All absolute minima of S are at = 5ij e N <= . Here a is a permutation 
of 0, 1, 2, . . . , N c — 1. There are no other minima. 

For very large (3/N c , the relative weights of distinct minima are given by the ratios of determinants of small 
fluctuations around the configurations. We parametrize the fluctuating degrees of freedom by the off diagonal entries 
of a hermitian matrix H which generates a conjugation around a diagonal saddle. 

SUa = i[U a , H]ij = iH l3 \<T^r - e 2 ^] (29) 
The diagonal terms of U don't contribute, and the second order action is 

6*S = 16/3 £ \H V |» sin 2 sin 2 ^A. (30) 

Since 



n i si »^ai| = n i si »^»] (31) 



for any permutation a, all minima should be visited with equal frequency for large f3/N c . 

One may ask whether the barriers between different minima are high and whether they become infinite in the large 
N c limit. The minima can be connected by sequences of transpositions. We do not know whether the barriers are 
infinite at infinite A^ c and whether in two dimensions, where the Feynman diagrams capture all the physics, one can 
ignore the many minima. 

From now on, we will use b — as the parameter held fixed as N c goes to infinity. In the planar limit, one 

approaches the continuum limit (and hopefully the 't Hooft solution) when b — > 00. Comparing to a continuum 
definition of the gauge coupling where the partition function is given by 

cx p{-iE T <}. ( 32 ) 

we get 

92 = W = 2^' (33) 
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To get a feeling for the numbers, consider two-dimensional QCD with massless quarks. The lowest scalar meson 
mass is approximately 0.967/6. In taking a sequence of increasing /3's and iV c 's with b 2 = 5, we find that this mass 
approaches 0.43 in lattice units. With a mass scale this large in lattice units, we can expect good agreement with the 
continuum theory only for small momenta. 

In two dimensions, quenching is not necessary Jl0| , [39| |. We will use a version of the model where there is no longer 
a need to perform the average over the angles. In this version, the eigenvalues of the matrices are still frozen to 
the set consisting of the iV c roots of unity. This adds a new "translational" symmetry to the model in which each 
Tjt is replaced by z M T M with z^ c — 1. This transformation preserves the eigenvalue sets. With N c chosen to be a 
prime number, the eigenvalues of will span the same roots of unity for any < k < N c — 1, and will simulate the 
situation on a torus of size N 2 . 

Let us summarize the symmetries of the bosonic action: 

S g {Ti,T 2 ) = S g (z 1 T 1 ,z 2 T 2 ) with |^| = 1 
S g (T x ,T 2 ) = S g {Tl,T 2 ) = S a (T x ,Tl) 

S g {T 1 ,T 2 )=S g {T 2 ,T 1 ) (34) 

These symmetries reflect invariances of the original toroidal lattice theory. 



B. Fermions and meson observables 

The fermionic momenta p M can be chosen from the set -j^-, j — 0,1, ... , N c — 1 for fermions with periodic boundary 

conditions and from the set j = |, §,..., N c — | for fermions with anti-periodic boundary conditions. Thus, the 
momenta flowing in the fermion lines are identical in type to those flowing in the gauge boson lines. With these 
momenta, the planar Feynman diagrams are similar to those of a system on a lattice with N 2 sites. 
The symmetries of (|34|) imply that the meson field propagators (^3j) obey the following relations: 

G(P, P 2 ) = G(-P 1 ,P 2 ) = G{P U -P 2 ) = G(-Pi,-P 2 ) 

G(P 2 ,P 1 ) = G{P 2 , -P x ) = G(-P 2 ,P) - G(-P 2 , -P) (35) 

These symmetries will be strictly enforced by explicitly adding the contributions of all gauge field configurations 
connected by the symmetry transformations of (|34|). 

G will be denoted by S for the scalar meson propagator and by P in the pseudoscalar case, 't Hooft's solution, 
together with subsequent work fl40| , provide explicit formulae for S and P in momentum space. Since they are 
composite, there are ultraviolet divergences, and the unrenormalized formulae are divergent. Since their source is in 
free field theory, these divergences can be regulated in many ways. We define the renormalized meson correlators 
Sr and Pr by subtracting at zero momentum for massive quarks and at an arbitrary momentum point for massless 
quarks. 

The continuum results we strive to reproduce using our reduced model are expressed in terms of the eigenvalues 
and eigenfunctions of 't Hooft's Hamiltonian H . It acts on square integrable functions defined on the interval [0, 1]. 
The boundary conditions are Dirichlet for massive quarks and von Neumann for massless quarks. For massive quarks, 
H is positive definite. For massless quarks, it has one eigenstate with zero eigenvalue, and all other eigenvalues are 
strictly positive. There are no degeneracies, and the eigenvalues of H provide the masses squared of the meson bound 
states. In the following discussion, all dimensionful quantities are converted to their dimensionless counterparts using 
the scale 

e 2 ^ (36) 
which has the units of mass squared. In these dimensionless variables, 

H<t> n =nl<l>n, n = 0,1,2,.... (37) 



The quark mass enters by the dimensionless parameter 



i PS) 



The squared momentum P 2 is replaced by the dimensionless variable Q 2 



p2 

Q 2 = (39) 
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The basic formulae for the meson propagators with massive quarks are: 



n>l odd 



(40) 



Asymptotically [i^ ~ 7r 2 n, and the residues approach a finite limit, so the subtracted sums converge. 
The residues r 2 are given in terms of the real, normalized wavefunctions <fi n (x): 



dx 



<t>n (x) 



(41) 



For massless quarks, one needs to take the limit 7^0. The relevant combination is Sr + Pr, which is the regulated 
version of (|27]). For m q = 0, the lowest meson mass is zero, and there is an infrared divergence at p = 0. This is 
handled by making the subtraction at nonzero Q 2 . More information about H and the evaluation of ([l0]) is collected 
in the appendix. 



C. Numerical results in the reduced model. 



We have focused on two cases. The first is 7 = 1 which, as explained in the Appendix, corresponds to an intermediate 
quark mass. The other case is massless quarks, where the chiral symmetry is spontaneously broken at weak c/ 2 -coupling 
in the planar limit. Other values of 7 are briefly discussed. 

Numerical results are obtained by working in the "unitary gauge" and generating gauge fields T\ that are related to 
Y by conjugation. Every conjugation belongs to an SU{2) subgroup of SU(N C ). A proposed conjugation is accepted 
or rejected by the Metropolis algorithm. The SU(2) elements are picked so that the acceptance rate is close to 0.5. 
One iteration is defined as a sequence of attempted updates for an entire set of jVc ^ c ~ 1 ) SU (2) subgroups of SU (N c ). 
The subgroups in the set are defined by placing the SU(2) matrix in the (n,n), (n,m), (m,n), (m,m) entries of the 
SU(N C ) T\ matrix. The indices obey 1 < n < m < N c . Fermionic measurements are made every 100 iterations, and 
1000 iterations are used for thermalization. 

Fermionic propagators are computed using (p4|). This is done by diagonalizing H\y exactly, constructing e, and 
then performing an exact diagonalization of 736. The propagator, g(pi,P2), at any mass, is obtained from the spectral 
decomposition of 736. This is done for all values of p\ and P2 and at a fixed value of T\. The "bare" scalar and 
pseudoscalar propagators are 

So(P) = ^ E Tr[g{p 1 + P 1 ,p 2 + P 2 )g( Pl ,p 2 ) 

C Pl,P2 

+ g(pi - Pi,P2 + p 2 )s(pi ; P2) 

+ g(Pl + P2,P2 + Pl)g{pi i P2) 

+ g(pi + P2,P2-Pi)g(pup 2 )] (42) 



Po(P) = ^5 E Tr [3(Pl +Pl,P2 + P2)g\pi,P2) 

C Pl,P2 

+ g(pi -P%,P2 + P2)g\pi,P2) 
+ g(pi + P2,P2 + P\)g\pi,P2) 

+ g(p 1 +P 2l P2-Pi)g\piiP2)\. (43) 

There are different vectors P — (Pi, P2) which correspond to the same value of P 2 = 4^^ sin 2 -j*-. By definition, the 
lattice meson propagators are a function of P 2 only. This function is obtained by averaging over all momenta with 
the same P 2 . After a final average over several values of T%, these propagators are the reduced model's approximation 
to©. 

We first compare numerical results to the continuum keeping 7 = 1 fixed. In order to study the approach to the 
continuum limit, we carried out simulations at five values of b 2 : 1, 2, 3, 4, 5. At each value of b 2 , we generated gauge 
field configurations at N c = 31,37, and 41 to see how the limit N c = 00 is approached. The lattice momenta VP 2 
range from to 2^/2 while \J Q 2 , reaches a maximum of 7 at b 2 = 1 and a maximum of 16 at b 2 = 5. In this range 
of momenta, good estimates of Pr(Q 2 ,1) and Sn(Q 2 7 l) are obtained by including only the first 100 terms in the 
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FIG. 6: Plot of S L (Q 2 ,1) at iV c = 41 for b 2 = 1,2,3,4,5. 
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sum over states. Since we have included the first 2000 eigenvalues (see Appendix), our "theoretical" data are very 
accurate. 

Let Sl{Q 2 , 1) and Pl(Q 2 , 1) denote the regularized, reduced model scalar and pseudoscalar propagators. As in the 
continuum, we subtract the propagators at zero momentum. Data at b 2 — 5 are shown in Figs. [}] and ^. At lower 
values of b, the N c dependence is even weaker. Even at b 2 = 5, we find that N c = 41 is large enough to provide 
accurate numbers for N c — oo. 

We study the approach to the continuum limit in Figs. || and g. Looking at different b 2 values, we see that the 
pseudoscalar propagator approaches its continuum limit more smoothly than the scalar propagator. For \J Q 2 < 4, 
the numbers seem to have converged to their b 2 — > oo limit, and there is reasonable agreement with the continuum 
result. On the other hand, the scalar propagator has not yet converged to its b 2 — oo limit in any region of Q 2 . 
It shows some overall agreement with continuum results, but it is less convincing than in the pseudoscalar case. In 
addition, the pseudoscalar propagator seems to approach its continuum limit monotonically, but the scalar propagator 
does not. 

The main difference between the scalar and pseudoscalar propagators is that the lowest mass in the pseudoscalar 
channel is smaller than the lowest mass in the scalar channel. In the pseudoscalar propagator, we have agreement 
below \[Q 2 = 4, and the lowest pseudoscalar mass is (io = 2.69713, which is below \[Q 2 = 4. However, the lowest 
scalar mass is fix — 4.16036, which is slightly above yQ^ = 4. Hence, all the masses contributing to the scalar 
propagator are above the region of Q where the pseudoscalar propagator converged well to its continuum limit. Thus, 
the difference we found between the scalar and pseudoscalar propagators is likely to be a finite b 2 effect. 

The difference between the two propagators can also be seen by extracting an effective mass. We set the momentum 
in one direction to zero and Fourier transform on the other component to produce a function of a conjugate, discrete 
"spatial" variable. Figure [?] is a semilog plot of the scalar and pseudoscalar propagators at sample parameter values 
b 2 = 5 and N c =41. We see that the pseudoscalar propagator is a better fit to a straight line. 

One can get an effective mass in either channel by forcing straight line fits. The resulting effective masses at N c = 41 
are shown in Fig. || as a function of b 2 . Although both masses appear to approach the correct continuum values, the 
finite b 2 effects on the scalar mass are much stronger. 

We now turn our attention to the behavior of the pseudoscalar propagator as a function of m q /e. We numerically 
computed this propagator as a function of \J~P 2 at several values of m q /e for b 2 = 5 and iV c = 37. An effective mass 
was extracted as before. This effective mass is plotted as a function of m q /e in Fig. ||. 

The empirical fit 

!^ = 2 + 0.51(^)- - 78 (44) 
m q e 

shown in the plot works quite well. 

With increasing m q /e, the ratio m 7r /m 9 approaches 2, as expected. However, the subleading behavior does not 
match the expression derived in the Appendix which has —1.33 rather than —0.78 in the exponent. On the other 
hand, m^/e is clearly not proportional to \Jm q je as would be appropriate for light quarks. We get 0.78 instead of 
0.5 in the exponent. The fit indicates that the masses we used are in an intermediate region. 

Finally we look at massless quarks. To compare the lattice data to theory, we need to pick a subtraction point. 
We used x/Q 2 = 3. For any finite N c , the scalar and pseudoscalar propagators are identical for massless quarks. As 
explained earlier, the common value of these reduced model propagators ought to be compared with the continuum 
Mr defined by: 

n>0 

The comparison for the subtracted meson propagator is shown in Fig. |l^ and is quite good. The deviations that 
now appear at small Q can be mostly attributed to finite N c effects which arc sizable now because of the massless 
pseudoscalar mode. 

IV. SUMMARY AND CONCLUSIONS 

The main objective of this paper was to introduce a new program whose aim is the numerical solution of planar 
QCD. Our first test was an application to two-dimensional planar QCD. Much more can be done on this test case; 
our results are just the very first steps. We hope to carry out similar preparatory studies in four dimensions. 



(45) 
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FIG. 7: Plot of the scalar and pseudoscalar propagators in real space at b = 5 and N c = 41. 
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FIG. 9: Plot of the effective mass extracted from the pseudoscalar propagator as a function of m q /e at iV c = 37. 




FIG. 10: Plot of the meson propagator for massless quarks as a function of b 2 . The b 2 = 1 data is at N c — 37, the b 2 = 2, 4 
data are at iV c = 41 and the b 2 = 3, 5 data are at N c = 47. 
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Our numerical experiment in two dimensions shows that reduction can reproduce some correct numbers, but that 
precise agreement would require a more substantial effort. It should be kept in mind that QCD 2 in the planar limit 
is quite challenging numerically. This model was chosen as our test case only because we have exact theoretical 
information about it. In four dimensions, the planar limit of QCD might be approached more rapidly. 

Perhaps surprisingly, two dimensions might be even harder than four dimensions in the reduced model framework. 
If all minima of the pure gauge action are correctly sampled, N c ~ V~3 . Therefore, roughly, the number of integration 
variables, N 2 decreases with the dimension at fixed equivalent number of lattice sites. One might not need more 
lattice sites in four dimensions than in two because infrared effects are stronger in two dimensions. 

The ultraviolet cutoff is increased by decreasing the 't Hooft scaled coupling e 2 , which is the same as taking our b 2 
parameter to infinity. Numerically we have learned that one needs larger N c values to see convergence to the planar 
limit for larger values of b 2 . The associated cost could be reduced by improving the lattice action and our fermionic 
operators. Already from [||, we learn that simple improvement methods have a beneficial effect. After reduction, the 
relative benefit from improvement might be less than in regular simulations. Of course, one needs to rethink what 
kind of algorithms are best suited for reduced models. 

We have determined numerically that the eigenvalue distribution of the parallel transporters around plaquettes 
has support in only a small neighborhood of unity when we are at reasonably weak gauge coupling. As a result, 
the spectrum of the Wilson Dirac operator develops a gap around zero. This in turn would speed up the numerical 
algorithms used to implement the lattice overlap propagators. 

Reduction would also provide a means for evaluating 1/N C corrections. Thus, one might not only get numbers 
in the planar limit but also an estimate for their accuracy in the context of real QCD. By going to the Veneziano 
limit instead of to the 't Hooft limit, we might get meson widths in addition to the meson masses. Recall that 
the momentum enters very differently here so that one might even contemplate a direct analytical continuation in 
momentum to physical values. The Veneziano limit, as a function of Nf/N c , would provide quantitative assessments 
of the valence approximation — something of great value at this point in time. 

Furthermore, the old problem of dealing with complex actions numerically might be more tractable in reduced 
models because the number of integration variables seems to be so much smaller. This would open new ways to look 
at ^-dependence and at finite baryonic density. 

We conclude that pursuing our proposal further is a risk that is worth taking and look forward to further develop- 
ments in this direction. 
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APPENDIX A 



In this appendix we describe the combination of numerical and analytical methods used to produce the theoretical 
predictions for the meson correlators. Without exception we are only dealing with the case of quarks of degenerate 
mass. Our presentation is fairly detailed. 

't Hooft's Hamiltonian, H, is defined by: 



(H<f>)(x)=j(- + -^—)<f>(x) 
x 1 — x 



dy 



P 



Mv) -<!>(?)] 



(Al) 



1. Asymptotics of wave functions 



There is an x — > 1 — x symmetry under which the eigenfunctions of H will be either even or odd. We wish to 
determine the correct boundary conditions on the wave functions and thus define the Hilbert space on which H acts. 
We shall focus on the x = endpoint. The above reflection about x = ^ determines then the boundary conditions at 
the other endpoint. 

We write = 
4>' {x) in the integral, ending up with 



in equation (Al) and perform an integration by parts. Then, we add and subtract 

y x 



1 x 1 — 



(1) - ct>{x) - 0(0) 



x 



2 \ x 



1 



II 



y-x 



We now assume asymptotics of the type 



<fi(x) ~ x 13 for x — > 



(A2) 



(A3) 



We plug this in and evaluate for x ~ 0. We use the asymptotic form over the entire integration range inside the 
integral, but shall later make sure that only the region of y close to zero contributes to the final expression. Also 
assume that < f3 < 1. Then the two most singular terms at x = are contained in: 



(H(f>){x) ~ (7 - l)^" 1 + ^-x - 1 \ogx 2 



-(3x 



0-i 



du 



/- 1 - 1 
it- 1 



u + 1 



+ Px 13 ' 1 log - 

X 



(A4) 



Integration variables were changed in the integral: y = ux 
integral is trivial and gives the second term that contains a logarithm in cq 



The term — ^ was added and subtracted 



(A4) 



The subtracted 
The integral converges as the 



upper limit goes to infinity. Thus, even if we stop the integration at a large finite number A the answer would be 
some number which will depend on A only slightly. But, if A is finite, the range of y's is from zero to Ax, so it is all 
in the asymptotic regime. This is why we were justified in using the asymptotic expression for <f) under the integral 
sign, for what seemed to be the entire range. The integral (for x = 0) is done using the n — ► 1 limit of 



r 

Jo 



du— 



7r(cot fiTT — COt /?7f) 



(A5) 



which is — 7r cot tt(3. 

The terms containing logarithms cancel and to cancel the subleading singularities (which are still too strong to be 
matched by a term linear in <p - which is what would be necessary for an eigenstate) we need: 



7 - 1 + 7T/3 COt 7T/3 = 



(A6) 



This formula can be found in 't Hooft's original paper 
For light quarks, 7 — > 0, and we get: 



0' 



37 for 7^0 



(A7) 
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For massless quarks the asymptotic analysis needs to be redone. In this case the Hamiltonian has the constant 
function as the lowest eigenstate. Thus, the action of H on a function is determined by the derivative of the latter. 

mix) = m^m - m ~ m - V(«)bg (—^ 2 

1 — x x 2 \ x 

y-x 
We assume 

<f>'(x) ~x a for x -> (A9) 

with < a < 1 and check for consistency. This will give a condition on a, which will have a solution in the required 
interval. As before, we take the asymptotic form throughout the integration range, and leave it to the end to check 
that the contributions kept only came from the asymptotic regime. We find: 

1 /"* ii a — 1 

(H(f>) (x) = x a + x a log x-x a du + • ■ • (A10) 

l + a J it - 1 

The u integral is rewritten by adding and subtracting terms chosen for their large u behavior, with the objective of 
isolating the contribution from the region near the endpoint y = 0. 



du = / du 

u — 1 



1 „_i 1 



u+1 



1 



log—- (All) 



Collecting terms in the integral and using 



duuP- 1 = - cot ^- for p < q (A12) 

1 — u q q q 



we get 



As a result, we get 



f 

Jo 



1 -2 

= 7rcot(7ra) (A13) 



1 x a 

(H<t>) (x) = + nx a cot ira + • • • (A14) 

a a + 1 



a is the unique solution between and 1 of the equation: 



7rcot7ra = (A15) 

a + 1 



This equation is the same as equation ( A6), extended to 1 < < 2, with [3 — a + 1 and with 7 = 0. Since the leading 
term is a constant, we see that the next period of the cotangent in ( |A6| ) provides the subleading behavior. 

In summary we conjecture that the eigenfunctions for massive quarks have the following expansion near x = 0: 

if>l{x) ~ {alx 131 + a\x^ +1 + •••) + {a 2 x fh + a\x^ +1 + • • •) 

+ • • • + (oj/' + a k ; X Pk+1 + •••) + ••• (A16) 



where k = 1, 2, • • • are, in increasing order, all the positive solutions of equation (A6). When the quark mass is 
taken to zero the entire series associated with the first root is replaced by a constant. 

In the massless case we also can say that, for any excited eigenstate of H, [H4>e){x) = E4>e(x), E > we have, 

4>e(x) = M ( 1 x a+1 ) + terms that vanish faster as x — > (A17) 

V a + 1 / 

The value of < a < 1 can be easily calculated numerically. Note that the subleading correction at the endpoint 
holds in an ever decreasing range as the energy of the state, E, increases. 

In conclusion, the Hilbert space appropriate to the Hamiltonian for the case of massless quarks consists of square 
integrable functions on the segment obeying von Neumann boundary conditions 4>'(0) = 4>'(1) = 0. In the massive case, 
the boundary conditions were Dirichlet, 0(0) = 0(1) = instead. Note that the internal Hilbert spaces corresponding 
to massless and massive quarks are distinct, as one might have guessed would be appropriate for a manifestation of 
spontaneous symmetry breaking in the planar limit. 
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2. Massive quarks: 7 = 1 



Equation (A6) simplifies when 7=1: its roots are 0k = k— i with fc = 1, 2, The double series in (A16) collapses 

into a single series and 4>n( x ) has the simple structure of \fx times a series in x for the x — » asymptotic regime. For 
the eigenvectors, there is a a; — ► 1 — ir symmetry, so the structure should be 



4> l n {x) = 2^x(l-x)f n (x) (A18) 

where /„ is either even or odd about the center of the segment (0, 1). It is convenient to map this segment to (—1, 1). 
For 7 = 1 the Hamiltonian is particularly simple: 

(Hct>)(x) = - [ dy—4f(y) (A19) 

Under the change of variables x — > 2x — 1, and still using a prime to denote derivatives with respect to the new 
variable we get: 

(Hcj>)(x) = -2 f dy—^(y) (A20) 

J-i y-x 

</> vanishes at the endpoints, so </>=constant is not an eigenstate. The eigenfunctions are defined over the segment 
(— 1, 1) and have the structure (dropping the superscript denoting the value of 7) 



4> n (y) = y/1 -y*f n (y) (A21) 

We shall use Chebyshev polynomials of the first (T) and second (U) kind to parametrize /„ jil]]. Some of their 
properties are listed below: 

U n (co S 0) = sisM£ 
T n (cos9) = cos n# 

f-i d y^V 1 -y 2u n(y) = -nT n+1 (x) 

J_ x dx^/l - x 2 U n {x)U m {x) = §£„, m 
2(x 2 - l)U m -x{x)U n - X {x) = T n+m (x) - T ]n _ m] {x) (A22) 

The indices n and m above are non-negative and vary in the ranges in which the equations make sense. 
We represent a given function <fi(y) as 



m = ^i-y 2 Y, a ^y) ( A23 ) 



n=0 

A short calculation then gives 



^ 00 

4>'{y) = —==Y,{n + l)a n T n+l {y) (A24) 

V 1 - V n=0 



Using the properties of the Chebyshev polynomials listed in (A22) we find 

(Hct>)(x) = 2 f dy— 7 =L= YVn + l) a „T n+1 (y) 

•'-l ^"V 1 -!/ n=0 

OO 

= 2?r^(n + l)a„Z7„(x) (A25) 

Suppose <p(x) is an eigenstate of H with eigenvalue £7, (H<fi)(x) — E<f>(x). Then, 

00 00 

E\/i -y 2 J2 a My) = 2tt J2( n + l)anU n {y) (A26) 



n=0 n=0 
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One can proceed now by using orthogonality on one of the sides of the equation. We choose to multiply the equation 
by \/l — y 2 U m (y) and integrate over y: 

E^a n dy(l-y 2 )U n (y)U m (y) = TT 2 (m + l)a m (A27) 

n n J —1 



We need to evaluate the integral on the left hand side. Using an identity in ( A22 ) we find 



(1 - y 2 )U n (y)U m (y) = - 
-l 1 



By y — > —y we see that £„ m = for n — m odd. Also, 



n=0 



\T\n-m\{y) — T n+m+ 2(y)] = £ nrn — 



— (m + l)a r , 



We defined new coefficients 



We also defined a new matrix X 



b n = vn+1 



1 



I —^/ILIL I — 

V Tn + 1 yn + 1 
We end up needing the eigenvalues of the symmetric matrix X: 



n=0 



^ ^ A mn 6 n — ~E^ 



To calculate the entries £ m „ we go to angular variables: y — cos 9: 

1 



d9 sin#[cos(n — m)6 — cos(n + m + 2)9] 



We need one elementary integral, for even integer k: 

d9 sin 9 cos k9 = — 

For odd k the integral vanishes. This leads, for n — m even, to 

1 



k 2 - 1 



In particular, on the diagonal we have 



(n + 771 + 2) 2 — 1 (n — m) 2 — 1 



1 



4(n+ l) 2 - 1 

For n — m odd the matrix element vanishes. The matrix X is given by: 



f 



y/(m + l)(n + 1) 



f 



(n + r7i + 2) 2 — 1 (n — m) 2 — I 



(A28) 



(A29) 



(A30) 



(A31) 



(A32) 



(A33) 



(A34) 



(A35) 



(A36) 



(A37) 



for 77 — 777 even and X mn = for 77 — 777 odd. If A is an eigenvalue of X, E — \- is an eigenvalue of the 't Hooft 
Hamiltonian. 

For high states one can use the diagonal terms as an approximation. One gets the asymptotic estimate 



E n ~7T 2 (77+1) 



(A38) 
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with the state label n starting at n = 0. 

We need to determine the normalization convention on the infinite vectors (£>0) &i> • ■ ■ > ) that would make the 
eigenstates 4>e(x) {H(f>E — E(/>e) normalized to unity by 4> 2 E {x) — 1. </>_e(x) is parametrized by 

00 uE 

4> E {x) = 2y/x(l - x) V —^=U n (2x - 1) (A39) 
where x is in the original range (0, 1). Introducing the expansion we find 

2 00 

i = ^E( & «) 2 ( A4 °) 

n=0 

where, 

oo 2 

^2x mn b n = ^b m (A41) 

n=0 

We are now ready to compute the residues for 7 = 1. We need to calculate 

p\ = \ -Ux) (A42) 



Here, 7 was set to unity and the eigenfunction is assumed correctly normalized. Changing variables, we arrive at 

^ = f>i n) 4 (A43) 

fc=0 

with 

h= d6 |sin(fc + l)6> (A44) 

Jo cos I 

We change the integration variable 9 to tt — 9 and average the expressions. Then, for k even we have to calculate 

h = / d9 — \ A45 
Jo sinS 

Use 

Mk+l)0 = et{k -2j)0 (A46) 

sin6> f-( v ; 

to conclude that for even fe only the j = | term in the sum makes a contribution, giving 

h = 7T (A47) 

For fc odd 

V^ + ^cosfl (A48) 

By the same technique as above we get non zero contributions only from j = (for odd k, k > 1). Hence, for k 
odd: 

h = -7T (A49) 

Only the absolute value of p x n is determined, since the sign of the wave function is not fixed by normalization. The 
final answer is: 



\p\\ = n 



k>0, k—n— even 



E -^= ( A5 °) 
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Here, the normalization condition is: 



E i^ f = ^ (A51) 



fc>0, k—n— even 

In each sector the eigenvalues were obtained numerically by diagonalizing a truncation of the infinite matrix X 



of (A31). T he re sidues were obtained using (A50) with the eigenfunctions of the truncated matrix X normalized 
according to ( A51). In each sector, a matrix of size 2000 x 2000 was diagonalized and this produced accurate estimates 
of the lowest 1000 eigenvalues in that sector. The lowest pole contributing to the pseudoscalar is at [iq — 2.69713 
and the lowest pole contributing to the scalar is at /xi = 4.16036. Asymptotically, the eigenvalues are given by 
/i^ = ir 2 n + jir 2 + O—. The residues, r^, reach an asymptotic value of n 2 . The even residues approach this number 
monotonically from above (tq = 11.7596864) and the odd residues approach this number monotonically from below 
(rf = 9.03541704). 

3. Massless quarks 

We now turn to the evaluation of eigenvalues and residues in the case of massless quarks. The correct boundary 
conditions for the massless case are von Neumann so we choose to diagonalize the Hamiltonian in the basis of cosine 
functions. 

(x\n) = s/2cosimx n > 1; (sc|0) = 1 (A52) 

The states are numbered, starting with 0. The ground state is known exactly and by keeping n > 1 we are in a space 
orthogonal to it. Just as before, even and odd n states do not mix, so we can diagonalize the Hamiltonian in each of 
these subspaces separately. 

(n\H\m) = f dx L dy / °° tdte~^ x ~ v ^ [cos(7rna;) — cos(7rny)] 

[cos(7rmx) — cos(7rmy)] (A53) 

To do the integrals we introduce, with no restriction on the integers n, m the integral 

I n , m =- f dx I dy I tdte-^ x - y \e mnx - e mnv \[e mmx - e mmv \ (A54) 
2 Jo Jo Jo 

By changing variables x — > 1 — x and y — > 1 — y we see that 

In,m = (— l) n+m I n ,m (A55) 

We now restrict ourselves to even n + m. This is enough to get all nonzero matrix elements (n\H\m). In this case we 
learn that I n ,m is real. 

Jo Jo 1 io°° tdte-^-^ [e mx - e mn ^ ] [cos(Trmx) - cos(Trmy)] (A56) 
But, since I n .m is real we can take the real part to prove 

(n\H\m) = I n<m + I n _ m (A57) 
We need only I n ^ m for n + m even. Using the symmetry under exchange of x and y we get 

[>1 f-X />oo 

In,m= dx dy tdte -t{x~y)^nx _ ^ny^rax _ ^ my ^ 



Using the above, and restricting in the end to n, m > 1 and even n + m we get: 

(n\H\m) = TT 2 nS n _ m + 2 tdt 

/ fl _(_ f_iy»-l„-tl r 1 I 1 ir 2 (n 2 +m 2 ) 

- (1 - e - *) J ^ n + ttt — tt ral ) (A59) 
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The matrix elements of H can be expressed in terms of sine-integral and cosine-integral functions: 

For n ^ m, with n + m even : 
(n|i/|m)=log(Bi^) 2 + ^±^log(f) 2 (A60) 



^i-rrfi [n 2 ci(irn) — m 2 ci(irm)] — 2{ci[ir(n + m)) + ci[rr(n — m)]} 

For n = m : 
(n\H\n) =7r 2 n-21ogn + log^-2(l+7) 



+2[imsi[im) + (-1)"] + 4ci(nn) - 2ci(2im) (A61) 



One need s to distinguish between Si, Ci and si, ci. Here are the definition for the functions used in equations (|A6 



and (A61) 



The asymptotics go as follows: 



si(x) = - /~ dt^ = -f + Jo dt B -if± = - f + Si{x) (A62) 
d( x ) = - dtsai = Ci(x) (A63) 

For integer and large positive k : 

s i( 7T k) = t^ + 0( J ^ s ) (A64) 
To find the residues in the case of massless quarks one needs to take the 7 — > limit on the massive formula: 

pl = Vl I' dx^- (A65) 
Jo x 

<\P n (x) is the n-th (ordered by eigenvalue and starting from n = 0) normalized eigenstate of the Hamiltonian with 
mass parameter 7. 

The entire contribution comes from the lower end of the integral, because only from there does one get a singularity 
in 7 that can compensate the vanishing prefactor. Thus, only the asymptotic behavior of the wave function is needed. 
It is given by: 

</>^(x) = A^x^^ + terms that vanish faster as x — > (A66) 

Hence, 



P° n = (A67) 

In general, is difficult to obtain, as it is fixed by the normalization of the wave function and depends on its values 
throughout the interval. But, for n = we know that, at small 7, cf>o(x) tends to 1. Hence, A® = 1 and we get 

Po = (A68) 

For arbitrary n we need to compute numerically the n-th wave function for the massless case, normalize it, and get 
A® t from its values at the end point. 

Al = 0°(Q) (A69) 

A sign ambiguity remains, as the sign of the wave function remains undetermined by the normalization condition. 
But, only (p^) 2 enters in the amplitude so everything is wel l det ermined. 



The eigenvalues were obtained from a truncation of H in ( A53 ) to s ize 1 000 x 1000 in each sector. The residues are 
obtained from the end point of the normalized eigenvector as given by ( A67). The lowest non-zero eigenvalue is a scalar 
meson (having an antisymmetric wave function) at [i\ = 2.4233 and its residue is r\ — 7.58941916. Asymptotically, 
the eigenvalues are described by /1 2 = 7r 2 n — 2 log(n) + 0(1). The constant term is close to ^7r 2 but due to the log(n) 
term it is hard to determine it accurately. The residues are relatively difficult to estimate due to slow convergence. 
They asymptote to r 2 = 7r 2 and are monotonic in n. 
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4. Light versus heavy quarks 

We wish to determine whether 7 = 1 should be viewed as a light quark mass or as a heavy quark mass. To do 
this we first find the leading behavior of the pion mass for very small quark masses. This behavior has the typical 
structure induced by spontaneous chiral symmetry breaking. We then compare the leading chiral approximation to 
the pion mass with the exact value. If the quarks are light the approximation should work well numerically. If the 
quarks are heavy the approximation should be off. We also compute the pion mass for very heavy quarks. Again we 
compare the approximate expression to the exact value. 

To compute the pion mass (n — state) to leading order in the assumed small quark mass we start from: 



7\2 



dxf {x){Hf ){x)= 1 [f {x) 



1 



tlx 





Symmetrize both terms, the first under parity x 



dy 



P 



1 -x 



(y - x) 2 

x, the second under interchange of x with y to get 



(WO 2 = 2 7 



<t>Z{x)¥ 



dx 



dy 



y-x 



(A70) 



(A71) 



As 7 — > the first term is dominated by the endpoint contribution. Actually, also the second term is dominated by 
the end point contribution. As before, we use the asymptotic behavior throughout the integration range, but ascertain 
at the end that we were dominated by the endpoints. From the first term we get a contribution: 



(A72) 



The second term requires more work: 



dx / dy 



!J" 



y-x 



dxx" 



du 



(A73) 



The u-integral is done extending the upper limit to infinity. The correction needed to account for this approximation 

r 1 2 

will be of higher order in (3. The u integral is done as follows: Write 
to obtain 



= — «T" u^T anc ^ ^° ^ ne integral by parts 



du 



,0 



1 



-i 2 



-1 - 2/3 



Jo 



,2/3-1 



,/5"l 



du- 



-1 - 2^7r[cot(2/37r) - cot(/37r)] = -(/3tt) 2 



0(/3 4 ) 



(A74) 



Doing the remaining x-integral we learn that the double integral gives a contribution equal to that of the first term. 
In summary we get, converting to physical units: 



This is the result found by 't Hooft Jl2) 
as a mass-square scale as 



2m q \g\J^+0(m 2 q ) 



(A75) 

Equation ( A75 ) can be rewritten in variables rendered unitless by using — — 



/4 



m-n- fA-n 

m q ~ V7 — 



1.90 

1 



(A76) 



For 7 = 1 we found numerically that ^ = 2.7 establishing that 7 = 1 corresponds to quite light quarks, but still the 

0(m 2 q ) corrections to eq. (A7E) make a substantial contribution. 

The pion is, by definition, the lowest energy eigenstate of 't Hooft's Hamiltonian in the sector of functions symmetric 
under x — > 1 — x. 
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For 7 large the first term in 't Hooft's Hamiltonian dominates, so the pion wave function will try to minimize its 
contribution by being concentrated around the point x = ^. Making the wave function too narrow incurs a price 
from the second term. This suggests a variational estimate for the leading and subleading order in ^ . The trial wave 
function is 



H-if 



(A77) 



4>(x) is normalized to unity up to corrections that are exponentially small in A. 
The expectation value of 't Hooft's Hamiltonian in the state <f> is: 



i i 



4?(x) + 



2io Jo 



dy 



<j>(y) - 4>{x) 



(A78) 



Contributions from the endpoints of the integral in the first term are exponentially small in A and will be neglected, 
under the assumption (to be later justified) that A diverges as 7 — > 00. Thus we can calculate the contribution of the 
first term by infinite range Gaussian integration, obtaining 



4 7 +^ + 0(A- 2 ) 
A 



The second term is 



dy 



g 2 ( x 2) — g 2 



x-y 



The integral (excluding the prefactor) has a finite A — > oo limit: 

/ 2 

/CO pC 
du / 
-00 J —1 



dv 



(A79) 



(A80) 



(A81) 



The expression (u — v) 2 in the denominator is represented by an integral over 

to{e -t(u-v) from 

zero to infinity. The 

u, v integrals of each one of the terms can be done by Gaussian integration, and calculating the appropriate 2x2 
determinants one finds that the above integral is given by: 



2tt 



Jo 



1 



2tt 



Hence 



Extremizing on A gives 



(H) 



47 + ^ + VttX 
A 



(A82) 



(A83) 



r 



(A84) 



Thus, indeed A is large for 7 large and the wave function is close to a delta function around x 
pion mass square we find: 



i. For the unitless 



For large 7 we expect therefore 



t<T, 



m a 



47+ (16tt)3-73 



2(1 + 0.697^) 



(A85) 



(A86) 



Using this formula blindly for 7=1 gives ^ ~ 3.4 instead of the true value of 2.7. We see that the true value 
is midway between the leading answer corresponding to light pions (1.9) and the leading plus subleading expression 
valid for heavy quarks. Thus, 7 = 1 corresponds to intermediate gauge coupling (in units of quark mass). 



